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Abstract. We have measured depolarized light scattering in liquid benzene 
over the whole accessible temperature range and over four decades in frequency. 
Between 40 and 180 GHz we find a susceptibility peak due to structural relaxation. 
This peak shows stretching and time-temperature scaling as known from a 
relaxation in glass-forming materials. A simple mode-coupling model provides 
consistent fits of the entire data set. We conclude that structural relaxation 
in simple liquids and a relaxation in glass-forming materials are physically the 
same. A deeper understanding of simple liquids is reached by applying concepts 
that were originally developed in the context of glass-transition research. 



PACS numbers: 61.20.Lc,61.25.Em,78.35.+c,64.70.Pf 

1. Motivation 

On short time scales, all liquids show solid-like elasticity. This has been impressively 
illustrated by Brillouin scattering of X-rays [Q: on a THz scale, sound propagates in 
water with almost the same speed as in ice, more than twice as fast as on the kHz 
or MHz scale of conventional ultrasonic measurements. Such a cross-over goes along 
with a decay of structural correlations; it is called relaxation, and more specifically a 
relaxation when it leads from solid-like to liquid-like response (correlations decaying 
to zero). 

When a liquid can be supercooled far enough, a relaxation becomes critically 
slow, so that the material ultimately becomes a glass. The dynamics of glass-forming 
liquids has been studied in great detail. Stretching and time-temperature scaling have 
been identified as generic properties of a relaxation. Additional scaling laws have 
been uncovered by a mode-coupling theory (MCT). Originally proposed as a theory 
of the glass transition, MCT is now generally recognized to offer a unified description 
of microscopic and relaxational motion at comparatively high temperatures where a 
relaxation occurs on a GHz scale. 

In experiments undertaken in both the supercooled and the normal liquid phase, 
a relaxation and mode-coupling dynamics are found to evolve continuously across the 
melting point of the concurrent crystalline phase: on a GHz to THz scale, the molecular 
dynamics seems to be insensitive to whether the liquid's state is thermodynamically 
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stable or not. This leads us to the hypothesis that the molecular dynamics in the 
normal liquid state will be insensitive to whether the liquid can be supercooled or 
not. A pioneering paper on water studies of metallic melts || and our own 
experiments on several molecular liquids suggest that indeed a relaxation and mode- 
coupling effects also occur in liquids that cannot be supercooled into a glass. 

For a more detailed test of our hypothesis, we now study one such liquid in 
depth. We choose benzene, which presents the following advantages: (i) the molecule 
is structurally very simple and highly symmetric; (ii) it is stiff on all relevant time or 
frequency scales ||; (iii) many often studied glass formers are structurally related to 
benzene; (iv) benzene is an excellent light scatterer. 

To investigate the dynamics of benzene over a wide dynamic range, we have used 
depolarized light scattering. A neutron scattering study is currently under way and 
will be published later; preliminary results support the conclusions we draw from 
light-scattering. 

2. Context and Theory 

2.1. a Relaxation 

In zeroth approximation, relaxation may be modelled by an exponential decay of 
correlations. This ansatz goes back to Maxwell's theory of viscoelasticity; it has been 
elaborated for dielectric response by Debye ||. The underlying physics is mean-field 
like: one considers thermal motion of an individual molecule upon which all the other 
molecules exert just a constant friction. 

However, in glass-forming liquids a relaxation is found to be stretched: it is spread 
much more than an exponential decay in time or a Lorentz line (in our context: a Debye 
resonance) in frequency. Popular fit functions assume a fractional time or frequency 
dependence, as in Kohlrausch's stretched exponential 

$kW cxexp[-(Vr)' 5 ], (1) 

or in the Cole-Davidson susceptibility 

XcdM oc [(1 + ^2 7 ^yT)^ ^l . (2) 

While the relaxation time r depend strongly on the temperature T, the exponents 
Pk or /3cd vary only weakly. This can be seen as a consequence of time-temperature 
scaling: in a good first approximation, a relaxation has the form 

$(t;T)~*(t/r(T)). (3) 

Towards high temperatures, this scaling law has an obvious limitation: a relaxation 
can never become faster than the temperature-independent microscopic modes. 
Indeed, on heating glass- forming materials towards the boiling point the curves log r 
vs. T become flatter and flatter, so that r approaches, but never reaches the intrinsic 
time scale of microscopic motion . Depolarized light scattering in molecular liquids 
could clearly resolve an a peak and confirm its scaling up to temperatures far in the 
normal liquid phase ||, |J . 

2.2. Relaxation in Simple Liquids 

Historically, relaxation in simple liquids (on a THz scale), and a relaxation in highly 
viscous liquids (originally measured mostly on Hz to MHz scales) have long been 
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seen as two genuinely different processes 10 . In the most simple monatomic liquids 
like argon or sodium, on which theory-oriented textbooks |ll], O, [l3) concentrate, 
characteristic relaxation times are of the order of 10~ 13 to 10 2 s, which is not 
much longer than the mean time between collisions |14| ], Under such circumstances 
relaxation is closely mingled with microscopic motion, and it is impossible to obtain 
isolated experimental information on relaxation alone. On the other hand, results 
of scattering experiments and molecular dynamics simulations cannot be understood 
without taking into account relaxation. Therefore, experimental data are usually fitted 
by theoretical expressions that contain memory kernels built upon an ad-hoc model of 
relaxation. 

Such fits, however, are rather insensitive to the functional form of the memory 
kernel, and therefore one seldom went beyond assuming simple exponential relaxation. 
In some cases, when fits were judged unsatisfactory, a sum of two exponentials was 
used fl5|| ; this approach, though admittedly arbitrary, has recently been revived in 
the analysis of X-ray Brillouin scattering on liquid metals |l6|, [ij], [ll| . 

In the investigation of glass-forming liquids this double exponential approach has 
long been overcome by formulae like (|l|) or (||). Today, after a relaxation has been 
observed across the melting point and up to a GHz to THz scale, it is no longer 
justified to consider relaxation in glass-forming materials and relaxation in simple 
liquids as two different physical processes. So we are led to suspect that the stretching 
and scaling properties of a relaxation hold in principle even in argon and sodium, 
although an experimental verification will be extremely difficult Jl9[ | . 

For the time being we prefer to investigate a molecular, non glass-forming 
liquid, benzene, which in a sense is intermediate between monatomics on the one 
side and molecular glass formers on the other side: the benzene molecule is small 
and highly symmetric, in contrast to glass-forming liquids which necessarily have a 
more complicated structure as to prevent crystallisation. Yet we will find structural 
relaxation in benzene to be slow enough 1 1 to allow for a direct, unambiguous 
observation of its stretching. 



2.3. Mode-Coupling Theory 

As mentioned in the beginning, mode-coupling theory (MCT) ||c], ^l| provides a 
unified description of a relaxation and low- frequency vibrations. It is a microscopic 
theory, formulated as function of wavenumber q and time i, and built upon the static 
structure factor S(q) and the density correlation function & q (t) — S(q,t)/ S(q). It 
starts with the formally exact equation of motion 

o = n- 2 $ q (t) + *,(t) + f dt' M q (t - t')$ q {t') . (4) 
Jo 

The memory kernel M q {t) contains fast and slow fluctuations, M q (t) ~ M' {t) + m q (t). 
The fast component, approximated as M' q (t) ~ "f q S(t), can be shown to be irrelevant 
for the long-time behaviour. 

The basic idea of MCT is now to expand the slow fluctuations m q (t) in 
polynomials of density fluctuations, and then to factorize all terms into pair 
correlations. In lowest order one obtains the bilinear functional 

m q mt)}= E ^iPkSpW **(*)• (S) 

p+k=q 
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The coupling coefficients V^pk can be derived from the static structure factor S(q); as 
S(q), they vary slowly with state variables like temperature T or pressure P. 

In this way the dynamics is completely determined by a closed set of integro- 
diffcrential equations. Depending on the numeric values of V^pk, the & q (t) either 
decay to zero or arrest at finite values. The border line T C (P) which separates these 
two cases has been called the ideal glass transition. 

For T < T c , the density correlations arrest at a finite Debye- Waller factor 
— > oo) = f q {T), as expected for a glass. On the liquid side, with T > T c , 
the & q (t) slow down on approaching a plateau f q {T), but ultimately they decay to 
& q (t — > oo) = in a process which is easily identified as a relaxation. Since the 
derivatives $ g (t) become negligible at long times, Eq. (Q) immediately reproduces the 
time-temperature scaling ([|), with the corollary that the line shape of $ 9 may vary 
with q. To first order, solutions of MCT equations are consistent with the Kohlrausch 
asymptote ffl). 

On cooling towards T c , a relaxation times are predicted to diverge with a 
fractional power law in T — T c . A comparison with measured relaxation times and 
viscosities [p2[ shows that such a divergence does not correctly describe the glass 



transition 23 . Instead, T c is found to describe a cross-over that is typically located 
15-20 % above the conventional glass transition temperature T g : while the density- 
density coupling of Eq. (||) becomes ineffective at T c , other transport mechanisms, not 
covered by the theory, remain active at lower temperatures. This interpretation of T c 
is supported by various other experimental indications of a cross-over pj, |25|, . 

In liquids to which MCT has been applied in the past, the cross-over occurs at 
shear viscosities of the order of 10 1 to 10 3 Poise |2^]. We note that this is much closer 
to the viscosity of, say, water at room temperature (about 1CP 2 Poise) than to the 
glass transition (which, according to wide-spread convention, occurs at 10 13 Poise). 
We note also that the dynamic predictions of MCT are expected to work best not in the 
immediate vicinity of T c , but at somewhat higher temperatures where the concurrence 
of low-temperature transport processes can be neglected. 

This suggests that MCT should be tested as a theory that describes the dynamics 
of glass-forming liquids at rather low viscosities, far above the glass transition, in a 
slightly supercooled or thermodynamically stable state. In the present work we want 
to show that the restriction to glass-forming liquids can be omitted altogether. 

2-4- Applying Mode Coupling to Real- Life Liquids 

Taken literally, Eqs. (@) and (§) assume a liquid composed of identical, spherical, and 
stiff particles: only in this case all interactions between the particles can be derived 



from S(q). Extending the theory to mixtures is straightforward j27, 28 , but including 
orientational and innermolecular degrees of freedom poses extreme difficulties: the 
notational and calculational efforts required by models so simple as a liquid made of 
linear molecules [^9[ |3Cj , or a dilute solution of linear molecules in spheres |2) are 
intimidating. Thus, a MCT of molecular liquids is presently not available. 

On the other hand, the time and temperature dependence of MCT solutions is 
insensitive to most of the structural information hidden in V'qpk- Taking into account 
orientational or innermolecular degrees of freedom may lead to new classes of solutions, 
but at least in some types of molecular liquids the fundamental mathematical structure 
of equations (^) and (|B|) will remain dominant. 

In such cases, MCT solutions can be characterized by quite few parameters. 
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Close to T c , the analytical expansions of $g(t) — f q depends in lowest order on just 
one nontrivial line shape parameter A. Complete time correlation functions can be 
generated by numeric solutions of very simple MCT models: 

In the minimal F i2 model just one correlator $(f) and two coupling 

coefficients in 

m(t) = v^(t) + v 2 Mt)} 2 (6) 

are sufficient to obtain relaxational stretching and the ideal glass transition. With 
just one more correlator, one can generate spectra with arbitrary a relaxation 

strengths /|: a bilinear memory kernel J34|, |35|, |36| 

m s (t) = u„$(t)$,(f) (7) 

couples 3> s (i) to <&(£), whereas does not depend on $ s (i). 

The so defined two-correlator F\2 model is a physically meaningful tool for fitting 
experimental data. As such it has already been successfully employed in several 



studies of glass-forming liquids |37], |38|, |39j, |40| ; the theoretical background is explained 



especially in Ref. J4fj[] . A detailed numeric study has confirmed the stability of such 
fits El. 



3. Light Scattering Measurements 

Benzene (T m = 279 K, T b = 353 K) was bought from Sigma Aldrich (99.9% puriss. 
p. a.), unpacked under inert gas and sealed into a Duran cuvette. To our surprise, 
the sample could be supercooled to 258 K where it remained liquid for several hours. 
Data were taken at seven temperatures between 258 and 352 K. 

Light scattering experiments were performed using a grating double monochro- 
mator U 1000 and a six-pass Sandercock-Fabry-Perot tandem interferometer. In order 
to achieve stable operation at maximum resolution, both instruments are placed in 
insulating housings with active temperature control. The optics around the interfer- 



ometer has been modified as described previously 42 43| , Depending on the spectral 
range, the interferometer is used in series with an interference filter of either 150 or 
1000 GHz bandwidth that suppresses higher-order transmission leaks of the tandem 
interferometer jl4|, Q below 3% or better. The filters are maintained at con- 
stant temperature. To account for any drift, the instrument function is redetermined 
periodically by automatic white-light scans. 

In the present experiment, the slits of the monochromator are set to 30 - 60 - 
60 - 30 /im, resulting in a resolution (fwhm) of 7.5 GHz; data are only used above 
200 GHz. The interferometer is operated with mirror spacings of 0.4, 2.8, and 16.3 mm, 
corresponding to free spectral ranges of 375, 54, and 9.2 GHz; some additional data 
were taken with 7.5 mm (20 GHz). 

On both spectrometers, a near-backscattering geometry (172°) is used to minimize 
scattering from transverse modes. In depolarized (HV) interferometer measurements, 
the usual leakage from the acoustic modes is seen; these lines are about a 100 times 
weaker than in polarized (VV) scattering, but still up to about 3 times stronger than 
the continuous HV spectrum. Subtracting separately measured VV spectra [E7I, we 
could completely remove the contamination from the HV data. 

Intensity calibration is always a problem in light scattering. Best results 
were obtained by matching all data to the middle (54 GHz) spectral range of the 
interferometer. In this range intensities are reproduced after a full temperature 
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Figure 1. Susceptibilities of benzene measured by depolarized light scattering 
at temperatures T =258, 268, 279, 293, 310, 330 and 352 K. In the low frequency 
wing, one finds the white-noise slope v , indicated by the dashed line. The 
nontrivial relaxational and microscopic dynamics inside the rectangular frame 
is shown on an enlarged scale in Fig. El 



cycle within about 3 %. The temperature-dependent intensity mismatch of other 
spectral ranges is higher and attains up to 20 %. Part of the problem may be due 
to distortions of the optical paths within the cryostat, which are particularly severe 
when the spectrometer is operated with small entrance opening. In the present study, 
intensities can also be estimated from a fit which is normalized by construction (see 
Figs. H below). 

Finally, the spectra are multiplied with the detailed-balance factor, 

I(v) = I{u) cxp(-his/2k B T) , (8) 

averaged over energy-gain and energy-loss side, and converted from intensity to 
susceptibility 

X!i{v)=l{v)/n{v) (9) 
with the symmetrized Bose factor 

n<y ^ ~ exp(hv/2k B T) - cxp(-hv/2k B T) ' ^ 

In other molecular liquids, comparison with neutron scattering E3, |4^, ^l], [52| 

has shown that depolarized light scattering yields at least qualitatively a good 
representation of the dynamic susceptibility, and therefore we will interpret Xisi 1 ') 
in very much the same way as a susceptibility from incoherent neutron scattering. 
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Figure 2. Enlargement of the intermediate frequency region of Fig. [lj Structural 
a relaxation leads to a peak between 40 and 180 GHz. The shoulder at about 
2 THz is associated with microscopic ballistic motion. The flat cross-over between 
these two peaks is a signature of mode-coupling dynamics: the data cannot be 
explained as a simple superposition of a relaxation and harmonic short-time 
motion. Jyolid li nns are fits with the mode-coupling two-correlator F\2 model 
(see Fig. |, Sects. 2.3 and 4.5 and Table m). 



4. Data and Analysis 

4-1- Susceptibilities on Logarithmic Scales 

Fig. [| shows susceptibilities from depolarized light scattering for seven temperatures 
between 258 and 352 K. In studies of glass-forming liquids, measuring susceptibilities 
over several decades and representing them on double logarithmic scales were decisive 
steps in detecting nontrivial, stretched relaxation In the case of benzene, the 

same procedure, on the same absolute frequency scale, is less rewarding: too much of 
Fig. [I] is filled by an uninformative v 1 white-noise wing. 

Therefore we show the nontrivial part of our data in Fig. || on an enlarged scale: 
the strongly temperature-dependent dynamics between 15 GHz and 3 THz. With 
increasing frequency, the x"(. v ) begin to deviate from the white- noise limit \" 1/1 > 
reaching a maximum between 40 and 180 GHz which we will ascribe to structural 
a relaxation. A comparatively flat region leads over to a shoulder at about 2 THz, 
above which the susceptibilities strongly decrease. Above 5 THz we find an extended 
gap; innermolecular excitations are only expected above 12 THz ||. 

The whole scenario is compatible with the high-temperature limit of what has 
been observed in many glass-forming systems. We note that the picture does not 
change up to the highest accessible temperatures: little below the boiling point, the a 
peak is still separated by almost a decade in frequency from the vibrational shoulder. 




Figure 3. The curves in the main figure are identical to the mode-coupling fits of 
Fig. pi except for being shown in absolute units, fulfilling the x" / v sum rule. The 
inset shows the amplitudes Amc(T) by which these curves had to be multiplied 
in order to fit the experimental data; these amplitudes essentially represent the 
Pockels coefficient by which light scattering couples to the microscopic dynamics. 



4-2. Absolute Intensities 

The temperature dependence of the scattering intensity is surprising: the height of 
the a peak increases by about 15 on heating from 258 to 293 K; then it falls back 
and reaches at 352 K about the starting level. The apparent anharmonicity in the 
high-temperature, high-frequency limit is unexpectedly pronounced, though a similar 
trend has been observed in several other liquids. 

To disentangle possible causes of these anomalies, we take advantage of mode- 
coupling fits. The physical meaning of these fits will be discussed later (Sect. |4.5p ; 
for the moment we take them just as a smooth parametrization of the measured data 
— with one specific advantage: Since the mode-coupling susceptibilities are obtained 
by Fourier transform of the derivative of a time correlation function, they obey the 
x"{v)jv sum rule by construction. Therefore we can take them as representing our 
light scattering data in absolute normalization. 

Resulting curves are shown in ^: they are strictly identical to the mode-coupling 
fits included in Fig. — except that the latter are multiplied by an amplitude 
Ayyc to adjust them to the arbitrary experimental intensity. Thanks to the intrinsic 
normalization, Fig. ^ shows a highly regular temperature dependence. In particular, 
we see no longer indications for a softening of the microscopic excitation spectrum 
at high temperatures: in the high-frequency wing, up to the boiling point all 
susceptibilities coincide, as expected for harmonic motion. The a peak height increases 
steadily with T; only between 1 and 2 THz little experimental imperfections are visible. 
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Figure 4. Master curve, constructed by rescaling the light scattering data of all 
seven temperatures (same symbols as in Fig. ^) to a common a peak. The peak is 
more stretched than any of the usual fit formulae is able to describe: a Lorentzian, 
motivated Maxwell-Debye viscoelastic theory, is completely inadequate (dotted 
line). Kohlrausch's stretched exponential (Eq. hi with /?k = 0.73, dashed 
line) and the Cole-Davidson function (Eq. j^, straight line) hold at least up to 
above the maximum of the peak. At higher frequencies, the extremly flat wing 
approximately follows a power law v~ b with b ~ 0.13 (dash-dotted line). 



The temperature dependence of Amc is shown in the inset of Fig. |3[ Up to 
293 K the Amc (T) scatter somehow, then they decrease systematically towards about 
75 % of the low-temperature average. A similar decrease of depolarized scattering 
intensity has been observed in many other liquids. However, lacking a means of 
absolute normalization, it was never clear whether this decrease reflected a property of 
the sample or of the scattering process. Unstabilities of the experimental setup added 
to the difficulty. Fig. || it appears now that the decrease of scattering intensities at 
high temperatures is not due to the sample dynamics; it rather appears that Amc(T) 
reveals a temperature variation of the Pockels coefficient that couples light scattering 
to the microscopic dynamics. 



4-3. a Relaxation 

For a quantitative analysis of a relaxation, we first test time-temperature 
superposition. Using the frequency-space representation of Eq. (||), and allowing for 
a temperature dependent amplitude, we rescale our data onto a master curve. Fig. |] 
shows that the line shape is independent of temperature up to at least five times the 
peak frequency. 

The a peak is obviously stretched, as can be seen by comparison to a Debye curve 
(dotted line) . The data are far better described by one of the empirical expressions (P 
or (g). The Fourier transform of the Kohlrausch stretched exponential [Eq. ([!])] fits 
the master curve up to about twice the peak frequency (dashed line) with a stretching 
exponent /3k — 0.73. 
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Figure 5. Mean relaxation times (r) as obtained from Cole-Davidson fits to 
the light-scattering susceptibility in the a relaxation regi on (V left scale). Solid 
lines are fits by the Vogel-Fulcher function (Eq. |l^, Sect. g-j ). The temperature 
dependence agrees in first order with that of the shear viscosity r\ (o, literature 
data |54), right scale). 



The Cole-Davidson function (g) with /?cd — 0.33 describes the master curve 
to even higher frequencies (solid line). Consequently, we use Cole-Davidson fits to 
determine the mean relaxation time 

(r) = J dt$(t)/$(0)=/JaDT, (11) 

which in turn is used in the iterative construction of the master curve. 

None of the empirical fit functions is able to fully describe the extremely fiat 
high-frequency wing of the a peak. For about one decade, this wing roughly follows 
a power law v~ h (dash-dotted line in Fig. The exponent is about b ~ 0.13, the 
precise value depending on the choice of the frequency range. Such a power law is 
remi nisce nt of MCT, though it should not be taken literal as an MCT asymptote (see 
Sect. U below). 

4- 4- Relaxation Times 

In a next step, we investigate the temperature dependence of the mean relaxation time 
(r), as determined from the Cole-Davidson fits [Eq. (p"T|)]. 

In Fig. H (t) is plotted as function of temperature and compared to the shear 
viscosity rj. In a good first approximation, (r) and rj show the same temperature 
dependence. This completes the demonstration that the susceptibility peak under 
study is indeed due to a relaxation, in the same sense as in any glass-forming liquid. 
Furthermore, the approximate proportionality (r) cx rj can be used to extend available 
viscosity data J54| by more than 15 K into the supercooled phase, and by 9 K towards 
the boiling point. 

Furthermore, Fig. || shows a Vogel-Fulcher fit 

(^cxcxpf-^M (12) 
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Figure 6. For a more detailed comparison of rj and (t), we eliminate 
their common first-order temperature dependence by plotting quotients. Shear 
viscosity data are the same as in Fig. 3; the relaxation times (t) are interpolated 
to the corresponding temperatures by means of the Vogel-Fulcher fit. The 
comparison (a) of (r) and v/T is motivated by the Stokes-Einstein relation. 
However, figure (b) shows that (t) agrees slightly better with r] than with n/T. 
Both plots are in arbitrary units. 



to the mean relaxation times, with Tq = 94.2 K and Eq = 650.6 K. We abstain from 
any physical interpretation, since several decades in (r) are needed for a meaningful 
verification of (|l^); we just employ the fit as a tool for a more detailed comparison of 
(t) and 77. 

This comparison is performed in Figs, ^a and where we divide either rj or 77/T 
by the Vogel-Fulcher estimate of (r). A proportionality (r) oc 77/T is suggested by the 
Stokes-Einstein relation D oc T/rj: when the diffusion constant D is determined from 
a time correlation function, 

$0) oc cxp(-TVi) , (13) 
one has (r) oc D~ x and thus (r) oc 77/T. Of course the stretched a relaxation in 



benzene is not correctly described by Eq. (13). Therefore the theoretical grounds for 
assuming (r) oc 77/T are rather weak. 

And indeed, Fig. ^ shows that (r) is not proportional to 77/T, nor to 77, but 
something in between. Such a temperature dependence has been reported at least 
once before: in a neutron spin-echo experiment on the high-temperature dynamics of 
glycerol p5| . We therefore conclude: the mean relaxation time observed by scattering 
shows the same temperature dependence as the shear viscosity — up to a prefactor of 
0(T) which at present no theory is able to predict. 

4-5. Mode-Coupling Fits 

We now extend our analysis beyond the a peak, considering the full experimental 
frequency scale up to some THz. The theoretical reference is given by mode-coupling 
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Figure 7. Coupling coefficients used in the schematic mode-coupling fits (Figs. ^ 
and jj| Sect. [l^L. Numeric values are also given in Table [l| The coefficients v\ 
and v<2 of Eq. (fj) control the intrinsic dynamics of the liquid, represented by 
the correlation function 3>(t); the slave correlator $s(iL which represents the 
experimental observable, couples to <t>(t) via v s [Eq. (M)]. In agreement with 
the spirit of mode-coupling theory, we find decreasing coefficients with increasing 
temperature. 



theory. MCT is perfectly compatible with the scaling properties of a relaxation 
obtained in the two preceeding subsections. Additionally, MCT predicts that a 
relaxation has a rather flat high-frequency wing which leads over to the microscopic 
molecular dynamics. This is just what we see in Figs. || [|. 

In glass-forming liquids, mode-coupling analysis usually concentrates on a scaling 
regime, designated as fast j3 relaxation, which is located between a peak and 
microscopic frequencies. When the a relaxation becomes sufficiently slow, the dynamic 
susceptibility passes through a minimum, and for frequencies around this minimum 
simple asymptotic power laws are predicted. In benzene we find a susceptibility 
minimum at the two lowest temperatures. This once again supports qualitative accord 
with an MCT scenario, though the (3 relaxation regime is not sufficiently developed 
to allow a scaling analysis. 

The approximate power law v~~ h (Sect. 4.3, Fig. ||) in the high-frequency wing 
of the a peak is of the low-frequency asymptote of fast (3 relaxation; however, the 
exponent b ~ 0.13 implies a line shape parameter A ~ 0.98 which, though formally 
allowed, is highly unlikely to represent the true asymptotic value, which whenever 
reliably determined has been found to fall into range of about 0.65-0.8. We therefore 
think that v~ h represents not more than a transient: somehow related to the scaling 
properties of MCT, but not representing an analytical asymptote of the [3 minimum. 

Therefore, we use numerical instead of asymptotic solutions of MCT. Specifically, 
we use the two-correlator F\2 model, introduced in Sect. 2.3, which is defined by 
the equation of motion (^) [with the set {Q q } replaced by the pair {$, <I> S }] and the 
memory kernels (^|), (^). While Q(t) models the intrinsic dynamics of the system, 
$53 (i) shall be interpreted as the correlation function observed by depolarized light 
scattering. 

The model contains seven parameters: two frequencies il, il s characterizing 
ballistic short-time motion, two damping coefficients 7, 7 S representing fast 
contributions to the memory kernel in Eq. (^), and three coupling coefficients vi, 
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Figure 8. In continuation of Fig. |?|, the remaining parameters of the mode- 
coupling fitSjare shown: (a) the microscopic frequencies Q (■) and Q B (□, kept 
fixed) of Eq. kl and (b) the damping coefficients 7 (♦) and 7 S (0) representing the 
rapidly decaying part M' q (t) of the memory kernel. Since we are using THz units, 
the figure shows strictly speaking Q/2n etc. 
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Table 1. Parameters of the two-correlator F12 model, as obtained from the 
least-squares fits shown in Fig. ti. Note that vi, V2, and v s are dimensionless, 
whereas Q, 7, and 7 S are given in GHz. The amplitude Amc i s m the arbitray 
units of our depolarized light scattering experiment. The eighth parameter, 
C s /2tt = 1000 GHz, was kept fixed. 



V2 and v s . These parameters are all expected to vary smoothly and monotonously 
with temperature. An eighth parameter, the amplitude Amc is not part of the model, 
but needed to adjust it to the arbitrary experimental intensity scale; these amplitudes 
have already been discussed above (Sect. 4J2, Fig. ||). 

The inner loop of the fit routine calculates 4>(i) and < & s (i) by iteratively solving 
Eq. (||) in the time domain [[56| ^7], |58|. Then & s (t) is converted into a susceptibility 
by blockwise Fourier transform, using the Filon method. The so obtained Xs{ v ) are 
fitted to the experimental data. 
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Figure 9. Mode-coupling coefficients v\ and V2 as in Fig. ^jin the phase diagram 
of the F12 model. The solid line indicates the ideal glass-transition singularity. On 
cooling liquid benzene through its entire range of existence, this phase boundary 
is only little approached. 

This procedure is performed independently for each of the seven measured 
temperatures. In an attempt to reduce the number of free parameters we find that 
the microsopic frequency of the slave correlator can be kept at a constant value 
fl s /2ir — 1000 GHz. All other parameters are found to show a weak, regular 
temperature dependence with only minor deviations from monotonicity, as can be 
seen in Figs. and ||. Values are also numerically given in Table 1. 

The time dependence of a relaxation is essentially given by the coupling 
coefficients v\ and v-i- In Fig. || the values obtained from our fits are shown as points 
in a phase diagram. For large values of v\ and V2, the F12 model becomes a glass. 
The phase boundary corresponding to the idealized liquid-glass transition is indicated 
in the figure. 

In benzene, the coupling coefficients fall clearly in the liquid phase; with 
decreasing temperature the glass-transition singularity is only little approached. This 
correlates with the fact that the measured susceptibilities show only a very first onset 
of a fast /3-relaxation minimum. For the same reason it would not be meaningful to 
use asymptotic scaling laws that are based on expansions in T — T c . From the available 
V\, V2, it is not even possible to extrapolate a hypothetical trajectory by which benzene 
would approach the glass transition if it could be further supercooled. Therefore it is 
impossible to indicate a meaningful value of the asymptotic line shape parameter A. 

5. Conclusion 

We used depolarized light scattering to measure the dynamic susceptibility of liquid 
benzene. Four spectral ranges of two spectrometers were combined to cover frequencies 
from 0.5 GHz to several THz. White noise prevails up to 10 GHz (Fig. [I]). Depending 
on temperature, a relaxational maximum is attained between 40 and 180 GHz. The 
high frequency wing of this maximum is extremly flat, and extends up to about 2 THz. 
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In the supercooled state, the susceptibility passes through a slight minimum around 
1 THz (Figs, g, ||). 

Such a broad relaxation pattern cannot be described by exponential memory 
functions that underly conventional theories of simple liquids. Instead, our results 
look very similar to what has been observed in many glass-forming liquids. This 
confirms our starting hypothesis, and provides the basis for our quantitative data 
analysis. 

As in glass-forming liquids, the relaxational a peak is stretched; it is even more 
stretched than the common fit formulae are able to describe. Timc-tcmperature scaling 
is obeyed with high precision and up to the boiling point, contradicting certain glass- 
transition theories which assume that a relaxation becomes Debye-like in the high 
temperature limit (Fig. 

Within the accessible temperature range, the mean relaxation time (r) of benzene 
varies by more than a factor of 4, and it is roughly proportional to the shear viscosty 
r\ (Fig. [5]). This accord is not improved by applying the Stokes-Einstein formula 
according to whic h (r ) should go with rj/T rather than 77 (Fig. ||). Implications are 
discussed in Sect. 4.4. 

Our observations are fully compatible with mode coupling theory. Originally, this 
theory attracted attention because of its ability to model a density-driven transition 
into a nonergodic state. Very soon, however, it became clear that this singularity 
does not describe glass formation. Instead, it is now generally recognized that MCT 
describes liquid dynamics at relatively low viscosities. In several studies of glass- 
forming liquids, fits were extended above the melting point of the concurrent crystalline 
phase, which was found to be irrelevant for the molecular dynamics under study. In 
our present work, we push this evolution one step further by applying MCT to a liquid 
that can hardly be supercooled (actually, benzene can be supercooled by nearly 20 K, 
which came out as quite a surprise). In our experiment, we cover the full range of 
existence, up to 1 K below the boiling point. 

In this temperature range, we can no longer apply the asymptotic expansions that 
are used in most MCT studies of glass-forming materials. Instead, we use numeric 
solutions of the full mode-coupling equations of motion. An elementary model with 
two correlators and three coupling coefficients is sufficient for a satisfactory fit to 
our full experimental data set. All parameters show a smooth, physically reasonable 
temperature dependence (Figs. Previous mode-coupling studies on glass-forming 

samples had mostly concentrated on the asymptotic predictions for fast [3 relaxation. 
A phase diagram makes clear that this scaling regime is not accessible in benzene 
(Fig. §). In such a situation numeric solutions of a minimal mode-coupling model 
provide the most adequate description of dynamic susceptibilities on the GHz to THz 
scale of structural relaxation and microscopic motion [p9[ . 
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